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Electromagnetic  signals  which  propagate  through  strongly  dis¬ 
turbed  regions  of  the  ionosphere  can  experience  scattering  which  can  cause 
appreciable  amplitude,  phase  and  angle-of-arrival  fluctuations.  This 
report  considers  the  performance  of  a  space  based  radar  (SBR)  that  must, 
operate  through  a  highly  disturbed  propagation  environment  such  as  might 
occur  during  a  barium  release  or  after  a  high  altitude  nuclear  detonation. 
A  brief  summary  of  the  propagation  channel  characteristics  is. given  in 
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terms  of  quantities  that  are  important  to  SBR  design  issues.  Previous  results 
for  coherent  detection  performance  are  summarized  and  new  results  are  given 
snowing  the  effect 'of  noncoherent  integration  on  target  detection  performance. 
Both  conerent  and  noncoherent  detection  performance  can  be  seriously  affected 
in  a  scintillation  environment  if  scintillation  is  not  considered  in  the  radar 
design. 
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SECTION  1 
INTRODUCTION 


Electron  density  irregularities  in  the  ionosphere  can-  produce 
random  variations  in  the  amplitude  and  phase  of  a  propagating  wave,  even 
at  frequencies  in  the  GHz  range  (Pope  and  Fritz,  1971;  Skinner  et  a!., 
1971;  Taur,  1976:  Fremouw  et  al.,  1978).  These  rapid  variations  in 
signal  amplitude,  phase,  and  angle-of -arrival  are  called  scintillations 
and  are  often  observed  over  satellite  communication  links  through  the 
ambient  ionosphere  at  VHF  and  UHF.  Strong  scintillation  is  occasionally 
observed  at  frequencies  as  high  as  L-band.  Since  even  small  fluctuations 
in  received  power  can  cause  degraded  system  performance,  the  effect  of 
scintillation  must  be  considered  in  the  design  of  a  space  based  radar 
(SEP)  which  must  operate  through  an  ionospheric  channel. 

Worst  case  or  Rayleigh  amplitude  scintillation  is  likely  tc 
occur  if  the  ionosphere  is  highly  disturbed,  as  fnr  example  by  high  alti¬ 
tude  nuclear  explosions  (Arendt  and  Soicher,  1964;  King  and  Fleming,  1980) 
or  by  chemical  releases  (Davis  et  al.,  1974;  Wolcott  et  al,.,  1978). 
Increased 'electron  concentrations  and  the  irregular  structure  of  the  ioni¬ 
zation  can  lead  to  intense  Rayleigh  signal  scintillation  at  frequencies  as 
high  as  the  7-8  GHz  SHF  band  (Knepp,  1977).  Consequently,  the  effects  of 
scintillation  are  important  to  any  UHF  through  SHF  radar  system  that  must 
operate  through  an  ionospheric  channel  and  that  may  have  to  operate  in 
highly  disturbed  environments. 
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In  an  earlier  work  (Dana  and  Knepp,  1983),  the  impact  of  a 
severely  disturbed  propagation  channel  on  the  coherent  detection  perform¬ 
ance  of  a  space  based  radar  was  investigated.  It  was  shown  that  disturbed 
propagation  environments  characterized  by  Rayleigh  fading  create  addi¬ 
tional  design  issues  that  do  not  occur  for  S3R  operation  iri  an  ambient, 
undisturbed  propagation  channel.  This  report  is  a  continuation  of  the 
earlier  work  that  considered  the  detection  performance  of  an  SBR  which 
Utilizes  coherent  processing  within  a  dwell.  Here,  detection  performance 
is  obtained  for  SBR's  that  perform  noncoherent  integration  of  the  ampli¬ 
tudes  from  multiple  dwells  that  form  a  radar  scan. 


SECTION  2 

RECEIVED  SIGNAL  DESCRIPTION 


In  this  report,  it  is  assumed  that  a  space  based  radar  (SBR) 
must  operate  through  a  severely  disturbed  signal  propagation  channel  that 
is  characterized  by  severe,  Rayleigh  fading  over  a  one-way  propagation 
path.  By  one-way  Rayleigh  fading  is  meant  that,  if  a  constant  signal  is 
transmitted  over  the  channel,  then  the  first  order  amplitude  statistics  of 
the  signal  received  at  the  channel  output  are  characterized  by  a  Rayleigh 
probability  density  function.  Rayleigh  statistics  are  worst  case  for 
propagation  of  electromagnetic  signals  over  a  wide  frequency  range  in  many 
different  kinds  of'  random  media  including  laser  propagation  in  turbulent 
air  (Fante,  1975),  and  VHF  propagation  through  the  ionosphere  (Fremcuw  et 
al.,  1978)  and  through  striations  composed  of  barium  ions  (Marshall, 

1982).  In  a  propagation  channel  disturbed  by  a  high  altitude  nuclear 
detonation,  large  portions  of  the  signal  environment  are  characterized  by 
Rayleigh  fading  (Wittwer,  1980;  Bogusch  et  al.,  1981). 

RECEIVED  SIGNAL  FIRST-ORDER  STATISTICS 

In  fading  conditions  it  is  convenient  to  write  the  total 
received  signal,  Sr,  as 

sr  *  S0[a/<o>]S  (1) 

The  received  signal  is  written  as  the  product  of  three  factors:  S0,  the 
mean  signal  that  would  he  received  from  a  point  target  (including  the 
effects  of  transmitter  power,  range,  frequency,  etc.);  S,  the  fractional 


change  in  tne  signal  due  to  variations  caused  by  the  propagation  channel ; 
and  o,  the  change  in  the  signal  caused  by  target  fluctuations.  The  factor 
S0  contains  the  mean  signal  level,  therefore  <S>  can  be  set  to  unity. 
Since  the  mean  signal  level  has  no  statistical  variation,  the  received 
signal  fluctuations  may  then  be  expressed  as  the  product  of  two  terms,  the 
first  representing  the  effect  of  the  target  variations  and  the  second 
representing  the  effect  of  signal  scintillation. 

Target  Statistics 

In  this  report,  it  is  assumed  that  the  target  cross  section 
fluctuations  can  be  described  by  either  Swerling  1  or  by  Swerling  2 
statistics.  In  both  cases  the  radar  cross  section  has  an  exponential 
probability  distribution  (or  equivalently,  a  Rayleigh  amplitude 
probability  distribution.)  Swerling  1  statistics  apply  to  the  case  that 
the  cross  section  remains  constant  over  an  entire  scan  (a  scan  is 
comprised  of  one  or  more  coherent  dwells  or  pulse-trains)  of  the  target. 
Swerling  2  statistics  apply  to  the  case  that  the  target  radar  cross 
section  varies  independently  from  dwel 1 -t'-dwel 1  but  is  constant  for  all 
the  pulses  that  comprise  a  coherent  dwell.  These  terms  describing  the 
transmitted  radar  signal  will  be  clarified  later.  In  either  case,  the 
probability  density  function  that  describes  the  target  cross  section 
variation  is  given  by  the  expression 

p(o)  =  _L_  exp(-«i/<a>)  (2) 

,  <o> 

where  ^a>  .  is  the  mean  value.  Equation  2  fully  describes  the  first-order 
statistics  of  the  target  cross  section.  To  completely  characterize  the 
first-order  statistics  of  the  signal  at, the  radar  receiver,  one  addition¬ 
ally  requires  the  statistics  of  the  signal  fluctuations  caused  by  propaga¬ 
tion  through  a  disturbed  channel. 


Propagation  Channel  Statistics 

In  the  case  of  satellite  communications  (SATCOM),  the  one-way 
propagation  path  is  often  the  only  case  of  interest.  Tn  the  case  of  an 
SBR,  there  are  three  possible  propagation  geometries  that  are  potentially 
important. 

One-Way  Fading.  Under  bi static  operation,  the  radar  transmitter 
and  receiver  are  located  separately  and  thus  there  may  be  two  independent 
propagation  paths,  one  from  the  transmitter  to  the  target  and  the  second 
from  the  target  to  the  receiver.  For  the  case  that  one  of  these  two  paths 
is  free  space  (undisturbed  by  fading)  and  the  other  is  severely  disturbed, 
the  received  signal  power  is  characterized  by  a  Rayleigh  probability 
density  function.  This  one-way  fading  case  is  identical  to  the 
characterization  of  a  severely  disturbed  SATCOM  channel.  For  this  case, 
the  probability  density  function  of  the  fluctuations  of  the  power  caused 
by  the  disturbed  propagation  channel  takes  the  form 

p  (S)  =  -L  exp{  -S/<S>}  (3) 

<S  >  . 

where  <S>  is  the  mean  value. 

Bj static  SBR  Geometry.  By  the  bistatic  SBR  geometry  here  is 
meant  the  case  that  the  transmitter  and  receiver. are  located  aboard 
separate  satellites  so  that  the  signal  travels  through  two  independent 
one-way  propagation  channels  both  of  which  are  severely  disturbed  by 
scintillation  effects.  In  this  case,  the  received  power  is  simply  the 
product  of  the  power  resulting  from  each  independent  one-way  propagation 
path.  The  probability  density  function  of  the  product  of  two  independent 
exponential  distributions  is  easily  calculated  as  (Papoulis,  1965) 
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(4) 


p  (S)  3  JL  K  (2/S/<S>) 

D  <S>  0 

where  K0  is  the  modified  Bessel  function  (Abramowitz  and  Stegun, 

1965).  <S>  is  again  the  mean  received  power. 

Monostatic  SBR  Geometry.  For  the  case  • f  monostatic  SBR  opera¬ 
tion,  the  transmitter  and  receive^  are  colocated  so  that  the  signal  propa 
gates  twice  over  the  same  path  passing  through  identical  irregularities. 
In  this  case,  the  received  voltage  is  proportional  to  the  square  of  the 
voltage  after  one-way  propagation.  The  received  power  is  similarly  pro¬ 
portional  to  the  square  of  the  received  power  after  one-way  propagation. 
The  probaoility  distribution  function  of  the  received  power  then  may  be 
obtained  from  Equation  3  using  the  transformation  S  (monostatic)  3  S 
(one-way)  with  the  result 

P_(S)  3  -===  exp{-^S<S>}  (5) 

m  /2S<S> 

where  <S>  is  the  mean  vaTue. 

Cumulative  Distribution  of  Channel  and  Target  Fluctuations 

The  effects  of  fluctuations  in  the  propagation  channel  and 
target  cross  section  fluctuations  can  be  combined  in  the  form  of  a  proba¬ 
bility  distribution  function  that  represents  the  product  of  two  independ¬ 
ent  random  variables.  The  cumulative  distribution  function  of  the,  re¬ 
ceived  power  is  given  for  the  three  cases  under  consideration  as 


(6a) 


Pi(Sp)  =  1 

W"1 

W  *  1  *  /  exp{-u-2Sr/(u2SQ)  }du  (6c) 

o 

where  the  subscripts  1,  b,  and  m  refer  to  the  one-way  propagation  path, 
the  bistatic  case,  and  the  monostatic  radar  case,  respectively.  K0  and 
are  the  modified  Bessel  functions  (Abramowitz  and  Stegun,  1965)  and 
S0  is  the  mean  received  power. 

In  the  absence  of  scintillation,  the  received  power  is  that  from 
the  Swerling  1  target  alone  and  the  cumulative  probability  distribution 
function  of  the  received  power  is  the  integral  of  Equation  2 

P  (s  )  »  1  -  exp{-$  /S0 }  (7) 

no  scint'  r'  r  '  ' 

Figure  1  is  a  plot  on  probability  paper  of  the  cumulative  proba¬ 
bility  distributions  for  the  four  different  situations  considered.  The 
curve  denoted  "no  fading,"  is  simply  a  plot  of  Equation  7  and  describes 
the  situation  of  a  Swerling  1  target  that  is  usually  considered  in  the 
design  phase  to  determine  the  radar  power  requi rements.  Since  this  curve 
also  applies  to  a,  Ray lefgtr  ampl  itude  distribution,  it  is  valid  for  the 
case  of  propagation  of  a  constant  signal  through  strong  irregularities  as 
would  occur  for  transionospheric  satellite  communications. 


-  2/Sr/S0  K  i(2/Sp/So) 

SO 

-  /  u  KQ(u)exp { -4S p/ (u 2Sq) }du  (6b) 
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Figure  1.  Cumulative  probability  distribution  function  of  the 
received  power  for  a  Swerllng  1  target  and  various 
propagation  geometries. 

It  is  apparent  that  the  effect  of  scintillation  is  to  increase  the  proba¬ 
bility  of  deep  lades  beyond  that  expected  for  a  Swerling  1  target  alone. 
For  example,  a  Swerling  1  target,  with  no  propagation  fading,  gives  a 
probability  of  10"1*  for  fades  of  40  dB  or  greater.  For  monostatic  opera¬ 
tion,  the  probability  of  a  fade  of  40  dB  or  greater  is  0.02.  It. will' be 
shown  that  this  Increased  chance  of  deep  fading  caused  by  the  disturbed 
transionospheric  channel  causes  degradation  of  target  detection  perform¬ 
ance  of  an  SBR. 

RECEIVED  SIGNAL  SECOND-ORDER  STATISTICS 

The  above  discussion  fully  describes  the  first-order  or  ampli¬ 
tude  statistics  to  be  expected  after  propagation  of  a  radar  signal  through 
a  severely  disturbed  channel.  Note  that  this  description  is  valid  only 
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for  the  worst-case  situation  where  Rayleigh  statistics  characterize  the 
one-way  propagation  channel.  Many  other  statistical  descriptions  are  pos¬ 
sible  for  the  received  signal  in  less  severely  disturbed  environments. 
However,  the  signal  description  here  leads  to  the  most  severe  effects  on 
SBR  performance. 

The  second-order  fading  statistics  are  specified  by  the  correla¬ 
tion  function  of  the  received  complex  voltage.  For  the  case  of  one-way 
propagation  of  an  initially  constant  amplitude  signal  through  a  severely 
disturbed  ionospheric  channel,  the  autocorrelation  function  of  the 
received  voltage  is  given  as  the  two-position,  two-frequency  mutual 
'oherence  function  (Knepp,  1983). 

Signal  Decorrelation  Tine 

For  a  K"4  in  situ  power  spectrum  of  three-dimensional  ionization 
i -regularities  between  outer  scale  Lo  and  inner  scale  the  decor¬ 
relation  time  is  (Knepp,  1983) 

T!  *  /l  L0/tn(L0/ti)o^vL  (8) 

where 

vL  *  the  velocity  of  the  line-of-sight  through 
the  center  of  the  ionized  layer 

o1  *  2(re*)2L0UNe2  rad2 

X  ®  RF  wavelength 

rg  *  classical  electron  radius  (2.82  *  I0"15m) 

t  »  thickness  of  ionized  layer 
ANg2  *  variance  of  electron  density  irregularities 


The  signal  one-way  decorrelation  time  is  an  inverse  measure  of  the 
fading  rate  or  fading  bandwidth.  Large  values  of  t:  correspond  to  slow 
fading  conditions  and  small  values  correspond  to  fast  or  rapid  fading. 


For  the  bistatic  SBR  propagation  geometry,  the  decorrelation 
time  of  the  received  signal  is  related  to  the  decorrelation  time  for  each 
of  the  individual  propagation  channels  according  to  the  relationship  (Dana 
and  Knepp,  1983) 


T 

0 


2  2 
T1T2 


1/2 


where  Ti  and  are  the  decorrelation  times  of  the  two  individual 
one-way  propagation  paths. 


For  the  monostatic  SBR  geometry,  the  two  channels  are  coincident 
(by  the  principal  of  reciprocity)  and  the  decorre’ation  time  is 

t0  *  x^/2  (10) 

where  ii  is  the  decorrelation  time  of  the  one-way  propagation  path. 

For  worst  case  Rayleigh  fading  the  correlation  function  of  the 
received  complex  voltage  always  has  the  Gaussian  form 

<E*(t+t)E(t )>  »  exp{-T2/t2}  (11)  ' 

for  any  q f  the  three  possible  geometries  (Dana  and  Knepp,  1983).  In  the 
following,  the  signal  decorrelation  time  will  always  be  denoted  by  the 
symbol  t0  independent  of  the  various  propagation  paths.  The  actual 
value  of  t0  is,  of  course,  a  function  of  the  geometry  as  well  as  the 
irregularity  structure  of  the  ionospheric  propagation  channel. 


In  this  report  it  is  assumed  that  t0  is  large  with  respect 
to  the  duration  of  the  transmitted  pulse  so  that  the  received  signal  is 
coherent  during  this  brief  time  period  that  is  typically  of  the  order  of 
several  tens  of  microseconds. 


Channel  Coherence  Bandwidth 


The  channel  coherence  bandwidth  is  a  measure  of  the  maximum 
bandwidth  available  in  the  propagation  channel  over  which  it  is  possible 
to  transmit  a  signal  without  imposition  pf  undesired  pulse  distortion. 

That  is,  in  a  fading  environment,  signal  spectral  components  separated  by 
less  than  the  coherence  bandwidth  exhibit  perfectly  correlated  fluctua¬ 
tions.  If  the  signal  spectral  components  are  separated  by  an  amount 
greater  than  the  coherence  bandwidth,  then  different  spectral  components 
will  undergo  uncorrelated  fading.  This  distortion  in  the  received  signal 
spectrum  causes  the  received  time  domain  signal  to  display  undesired  time 
sidelobes.  On  the  other  hand,  a  propagating  pulse  remains  undistorted  as 
long. as  the  maximum  instantaneous  signal  bandwidth  is  less  than  the  coher¬ 
ence  bandwidth.  This  is  often  the  case  for  frequency  hopped  radar  sys¬ 
tems.  However,  frequency  hopped  signals  that  are  separated  in  frequency 
by  an  amount  exceeding  the  coherence  bandwidth  will  experience  independent 
fading. 

For  a  K-4  in  situ  power  spectrum  of  three-dimensional  ionization 
irregularities  the  coherence  bandwidth  is  givpn  by 


/2  w2c(z  +z  )l 
.  t  r  o 

coh  r2>4tn(L  /e  )zz  LanJ 
e  o  i  t  r  e 


(1?) 


for  the  case  of  a  monostatic  SBR  propagation  geometry.  This  value  is 
smaller  than  the  coherence  bandwidth  for  one-way  propagation  by  a  factor 
of  1//2  (Knepp,  1982).  In  this  expression 


15 


c  =  velocity  of  light  in  vacuum 

zt  =  distance  from  the  transmitter  to  the 
center  of  the  ionized  layer 

zr  =  distance  from  the  target  to  the 
center  of  the  *oni zed  layer 

zt  +  zp  =  total  one-way  propagation  distance 

For  the  case  of  the  Mstatic  SRR  propagation  channel,  the  coher¬ 
ence  bandwidth  is  related  to  the  coherence  bandwidth  of  each  of  the  one¬ 
way  propagation  paths  by 

“cohl  “coh? 

“coh  *  —  ; 

“cohl  +  “coh2 

In  this  report.,  it  is  assumed  that  the  channel  coherence  band¬ 
width  is  sufficiently  large  compared  to  the  radar  pulse  bandwidth  so  that 
no  time  domain  distortion  of  the  received  pulse  occurs.  However,  varia¬ 
tions  in  the  coherence  bandwidth  with  respect  to  the  frequency  hopping 
bandwidth  of  a  jam  resistant  frequency  agile  radar  are  considered  here. 

It  will  be  shown  that,  such  variations  affect  the  probability  of  target 
detection  when  dwells  transmitted  at  different  radar  frequencies  are  non- 
coherently.  combined. 
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SECTION  3 

RADAR  SYSTEM  CHARACTERISTICS 


In  this  report,  it  is  assumed  that  an  SBR  operates  through  a 
disturbed  ionospheric  channel  to  detect  and  track  targets  near  the  earth's 
surface.  Thus,  relative  to  a  ground  based  radar  with  similar  functions, 
ar,  SBR  has  several  limitations.  First,  targets  are  detected  and  tracked 
at  very  long  ranges  Second,  available  onboard  transmitter  power  is  rela¬ 
tively  low. 


These  two  points  imply  low  received  signal-to-noise  ratio  per 
pulse  and  therefore  require  long  integration  times.  However,  the  cross 
section  of  a  moving  target  remains  constant  or  coherent  for  only  a  few 
tens  of  milliseconds  because  of  target  me  ion  and  resulting  constructive 
and  destructive  interference  between  many  scattering  centers.’  Hence, 
during  a  radar  scan,  the  total  energy  transmitted  at  a  target  is  divided 
into  a  number  of  dwells  (also  referred  to  as  bursts).  Each  dwell  consists 
of  some  number,  n,,cf  pulses  transmitted  at  some  radar  frequency,  which 
are  coherently  integrated  upon  reception.  The  radar  frequency  is  changed 
between  dwells  as  a  form  of  jamming  protection  and  to  assure  independent 
samples  of  the  target  cross  section.  The  detected  amplitudes  of  all  the 
dwells  which  form  the  total  radar  target  scan  are  then  noncoherently  com¬ 
bined  in  a  postdetection  integration  process. 

Figure  2.  shows  a  simplified  block  diagram  of  a  generic  SBR 
receiver.  The  complex  input  signal  contains  amplitude,  phase  and  doppler 
Information  from  target,  clutter,  and  thermal  noise  sources.  The  weighted 
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Figure  2.  Block  diagram  of  generic  SBR  receiver. 


discrete  Fourier  transform  (DFT)  simultaneously  acts  as  a  coherent  inte¬ 
grator  and  as  a  clutter  rejection  filter.  After  the  amplitude  from  a 
single  dwell  is  obtained  from  the  amplitude  detector,  the  radar  frequency 
is  changed  and  several  more  dwells  transmitted  at  different  frequencies. 
Any  number  of  these  independent  samples  may  then  be  noncoherently  combined 
and  the  resulting  amplitude  compared  to  a  threshold  and  a  "hit"  or  "miss" 
declared.  For  the  sake  of  simplicity  in  subsequent  calculations,  it  is 
assumed  that  the  target  doppler  is  zero  and  that  no  clutter  rejection 
filtering  is  applied. 
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SECTION  4 

PROBABILITY  OF  DETECTION 


In  this  report,  SBR  target  detection  performance  is  summarized 
for  the  case  of  one  dwell  per  scan  for  slow  and  fast  fading.  These 
results  are  described  in  detail  in  a  previous  report  (Dana  and  Knepp, 
1983).  The  case  of  multiple  dwells  per  scan  is  then  considered  in  detail 
and  results  are  presented  for  the  situation  in  which  a  frequency  selective 
propagation  environment  acts  to  decorrelate  the  propagation  disturbance 
from  dwell -to-dwell. 

PROPAGATION  EFFECTS  ON  COHERENT  INTEGRATION 

Since  an  SBR  will  be  required  to  operate  over  a  wide  range  of 
fading  rates  relative  to  the  coherent  processing  time,  a  distinction  can 
be  made  between  slow  and  fast  fading  that  is  an  invaluable  aid  in  under¬ 
standing  the  effects  of  scintillation  on  receiver  performance. 

Slow  Fading 

In  slow  fading  conditions,  the  duration  of  signal  fluctuations 
Is  very  long1  compared  to  the  coherent  processing  (or  dwell)  time.  Hence 
the  signal  amplitude  and  phase  are  relatively  constant  over  the  receiver 
coherent  integration  time  which,  in  the  expanded  receiver  model  of  Figure 
3,  is  the  time  to  Integrate  n  pulses.  In  this  case  only  dwel 1 -to-dwel 1  or 
scan-to-scan  signal  amplitude  fluctuations  affect  SBR  target  detection 
performance.  Therefore,  the  effect  of  slow  fading  on  target  detection 
performance  may  be  determined  solely  on  the  basis  of  the  first-order 
signal  amplitude  statistics  expressed  by  Equations  3-5. 


..  t 
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Figure  3.  Expanded  block  diagram  of  S8R  receiver. 


In  alow  fading,  the  amplitude  output  from  the  coherent  integra¬ 
tion  shown  in  Figure  3,  w  =  (u2+v2)1/2,  has  the  well  known  Rician  proba¬ 
bility  density  function  (Dana  and  Knepp,  1983) 

p(w|a)  a exp{-  {w2+n2a2)/2noN2}  I0(wa/aN2)  (14) 

n°N 

where  a  is  the  target  amplitude  per  pulse,  n  is  the  number  of  coher¬ 
ently  integrated  pulses,  2 o2  is  the  total  noise  power  per  pulse,  (o2  is 
the  noise  power  per  pulse  in  each  of  the  quadrature  channels),  and  I0  is 
the  modified  Bessel  function.  Equation  14  assumes  that  the  amplitude  a 
does  not  change  during  the  coherent  dwell. 

Now  consider  the  effect  of  target  cross  section  fluctuations 
while,  for  the  moment,  allowing  no  fluctuations  due  to  the  turbulent 
propagation  environment.  Thus  assume  that  S,  the  propagation  contribution 
to  the  power  as  given  by  Equation  1,  is  constant.  If  the  target  cross 
section  has  a  Swerling  1  distribution,  as  given  by  Equation  2,  the  target 
amplitude  probability  density  function  is  then 

p(ajS)  ■  2a  exp {-  a2/SS0}/SS0  (15) 

foi  a  fixed  value  of  S.  Now  the  probability  density  function  for  the 
envelope  w  may  be  calculated  as 
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(16) 


p(wlS)  =  j  p(w|a)p(a|S)da 
1  o'  1 

or  simply 

p(w|S)  -  - - - - - —  exP  |  5 - - - J  (17) 

1  noN2(l+n5S0/2oN2)  (  2n aN 2 ( 1+nSS 0/2aN2 ) I 

In  the  single  dwell  per  scan  case,  it  is  sufficient  to  consider 
only  the  amplitude  w.  After  detection,  the  signal  amplitude  is  compared 
to  a  threshold  value  and, a  hit  or  miss  is  declared.  The  threshold  value 
t  is  set  on  the  basis  of  the  noise  alone  to  achieve  a  desired  probability 
of  a  false  alarm  Pfa  given  by 

Pfa  ~  ^  p(wjs0=O)dw  =  exp( -t2/2noN2}  (18) 

t 

The  probability  that  the  signal  will  be  detected  for  a  given 
value  of  S  is  the  same  as  the  probability  that  the  signal  amplitude  w 
will  exceed  the  threshold.  Therefore 

1 , 

1 1 + < 3 NR  >S  ‘ 

Pd(S)  =  /  p(w Js )dw  =  PJ  ■  (19) 

where  <SNR>»nSo/2oN?  is  the  mean  signal -to-nqise  ratio  per  dwell.  Equa¬ 
tion  19  gives  the  well  known  relationship  between  false  alarm  probability 
and  probability  of  detection  for  a  Swerling  1  target  with  mean  power  $o 
for  a  fixed  propagation  condition  or  a  fixed  value  of  S. 

To  obtain  the  probability  of  detection  for  the  case  of  a  Swerl¬ 
ing  1  target  combined  with  propagation  through  strongly  turbulent  ioniza¬ 
tion,  it  is  convenient  to  calculate  th;  probability  of  detection  according 
to  Equation  19  for  a  given  value  of  power  S  and  then  to  average  over  the 
distribution  of  power.  The  probabi 1 ity  of  detection  is  then  written  as 
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1  , 

1+<SNR>S 


(20) 


P  -  /  P- 
d  1  fa 
o 


p  (S)dS 
Kpropv 


where  Ppr0p(S)  is  the  power  distribution  after  single  or  multiple  trans- 
ionospheric  propagation  and  is  given  by  Equation  3,  4,  or  5,  whichever  is 
appropriate  to  the  actual  geometry,  with  <S>  set  equal  to  unity. 

The  average  probability  of  detection  for  a  Swerling  1  target 
with  no  scintillation  and  with  one-way,  bistatic  and  monostatic  scintilla¬ 
tion  geometries  is  shown  in  Figure  4  for  a  false  alarm  probability  Pfa 
of  10-6,  typical  of  a  modern  search  radar.  For  the  one-way,  bistatic,  and 
monostatic  scintillation  geometries  the  results  are  obtained  using  numeri¬ 
cal  integration  techniques.  For  the  Swerling  1  target  with  no  scintilla¬ 
tion,  Equation  19  is  plotted  with  S=l.  As  was  to  be  expected  from  the 
cumulative  probability  distribution  statistics,  the  monostatic  case  has 
the  greatest  detection  sensitivity  loss  and  the  bi static  case  has  somewhat 
less  loss.  Relative  to  the  optimum  detection  curve  with  no  propagation 


Figure  4.  Probability  of  detection  as  a  function  of  SNR  for  slow 
fading  conditions. 


fading,  at  a  0.7  probability  of  detection  value,  the  detection  sensitivity 
losses  are  4  dB  with  one-way  fading,  7  dB  with  bistatic  fading,  and  11  dB 
with  monostatic  fading.  Thus,  it  appears  that,  during  severe  slow  fading, 
a  gain  of  4  dB  is  possible  by  using  bistatic  operation  instead  of  mono¬ 
static  operation.  This  4  dB  gain  in  detection  sensitivity  may  not  be 
realizable  when  the  differences  in  monostatic  and  bistatic  geometry  and 
radar  cross  section  are  taken  into  account. 

Fast  Fading 

In  fast  fading  conditions,  the  signal  decorrelation  time1  is  less 
than  the  S8R  coherent  integration  time  and  the  effectiveness  of  coherent 
processing  is  reduced.  Equivalently,  the  signal  bandwidth  is  spread  be¬ 
yond  the  receiver  bandwidth  so  that  signal  energy  is  lost.  This  loss  is 
caused  by  destructive  interference  of  radar  pulses  which  are  uncorrelated 
in  amplitude  and  phase  with  preceding  and  following  pulses  during  the 
coherent  integration  time.  This  fast  fading  imposes  additional  loss  in 
detection  sensitivity  beyond  that  imposed  in  slow  scintillation. 

Under  fast  fading  conditions,  the  analytic  approach  used  above 
is  no  longer  generally  sufficient  since  the  received  signal  is  decor- 
related  during  the  duration  of  a  dwell.  For  the  bistatic  and  monostatic 
SBR  geometry,  it  is  necessary  to  generate  Monte  Carlo  realizations  of  the 
received  voltage,  which  possess  the  required  first-  and  second-order 
statistical  description. 

To  compute  the  probability  of  detection,  many  independent  reali¬ 
zations  of  the  received  quadrature  voltages  are  generated,  each  realiza¬ 
tion  is  then  numerically  integrated  and  a  pseudo-random  sample  of  the 
integrated  noise  voltage  is  added.  At  this  point  the  integrated  voltages 
are  amplitude  detected.  This  signal -plus-noise  output  amplitude  is  then 
compared  to  the  threshold  and  a  hit  or  a  miss  is  declared.  The  detection 
process  is  repeated  many  times  to  obtain  probability  of  detection  statis¬ 
tics  as  a  function  of  signal  decorrelation  time. 
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Figure  5  shows  probability  of  detection  versus  the  mean  signal- 
to-noise  ratio  nSo/2aN2  in  dB  per  dwell,  for  monostatic  SBR  operation 
with  a  probability  of  false  alarm  of  10“6 .  In  the  figure,  the  probability 
of  detection  is  shown  for  a  Swerling  I  target  with  no  fading,  and  for  a 
Swerling  1  target  with  slow  fading.  These  curves  are  taken  directly  from 
Figure  4.  The  simulation  results  for  fast  fading  are  shown  as  points 
which  are  connected  by  smooth  dashed  curves  and  denoted  by  values  of 
-r o / Tc i ,  the  ratio  of  decorrelation  time  to  coherent  integration  time. 

To  obtain  these'  results,  each  dwell  is  formed  from  400  pulses.  .The  95 
percent  confidence  interval  is  shown  about  one  of  the  simulation  results. 


MEAN  SNR  (dB/ DWELL) 


Figure  5.  Probability  of  detection  for  monostatic  space  based  radar. 

Figure  5  shows  that  fast  fading  causes  .  .durable  additional 
loss  in  detection  sensitivity.  For  a  value  of  about  0.02  for  the  ratio  of 
signal  decorrelation  time  to  SBR  coherent  integration  time,  the  additional 
loss  is  about  8  dB  relative  to  the  slow  fading  Cro+*»)  limit  at  a 
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value  of  0.7.  As  the  decorrelation  time  increases  to  approach  the  coher¬ 
ent  integration  time,  the  values  approach  the  slow  fading  limit  as 
expected.  The  results  for  fast  fading  in  this  figure  show  that  the  slow 
fading  limit  for  P<j  is  exceeded  by  the  curves  for  values  of  T0/TCj  of 
0.09  and  0.18  for  mean  signal -to-noise  ratios  above  30  dB.  This  behavior 
is  the  result  of  two  competing  effects.  As  t0/Tc^  decreases  the  effi¬ 
ciency  of  the  coherent  integration  process  is  degraded.  However,  as 
T0/Tci  decreases,  the  number  of  independent  samples  of  the  received 
signal  increases,  thereby  changing  the  statistics  of  the  output  of  the 
coherent  integrator.  For  monostatic  operation,  with  the  very  severe  input 
amplitude  statistics,  the  change  towards  more  favorable  statistics  gives  a 
very  slight  gain  in  detection  performance  over  a  small  range  of  to/Tc-j 
values. 

PROPAGATION  EFFECTS  ON  NONCOHERENT  INTEGRATION 

In  this  section,  the  effects  of  fading  are  determined  as  it 
impacts  the  detection  performance  of  an  SBR  that  utilizes  noncoherent 
integration  of  many  dwells  that  comprise  a  scan.  Slow  fading  conditions 
are  assumed  in  the  remainder  of  this  report  so  that  the  signal  decorrela¬ 
tion  time  (due. to  fading)  is  large  with  respect  to  the  duration  of  all  the 
dwells  of  an  entire  scan.  Thus,  there  is  no  loss  due  to  pulse-to-pul se 
decorrelation  during  the  coherent  integration  process  shown  in  Figure  3. 

Furthermore,  this  slow  fading  assumption  introduces  the  feature 
that  the  propagation  effects  are  controlled  by  the  relationship  between 
the  radar  hopping  bandwidth  and  the  channel  coherence  bandwidth.  For 
example,  if  the  coherence  bandwidth  is  much  larger  than  the  hopping  band¬ 
width,  changes  in  the  radar  frequency  from  dwell -to-dwell  introduce  no 
changes  in  the  observed  scintillation  during  that  scan.  Conversely,  if 
the  radar  hopping  bandwidth  is  large  with  respect  to  the  channel  band¬ 
width,  the  voltage  component  due  to  scintillation  changes  from  dwell-to- 
dwell  providing  statistical  variation  and  giving  improved  detection 
performance. 
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First,  a  comparison  of  detection  probability  is  given  for 
ambient  environments  with  no  fading  for  the  case  of  Swerling  1  and  2 
targets.  Probability  of  detection  is  calculated  for  the  combination  of 
the  coherent  integration  of  the  voltage  from  n  pulses  during  a  dwell  and 
the  noncoherent  integration  of  the  amplitudes  from  m  dwells  during  a 
scan.  From  these  results,  it  will  be  seen  that  there  is  an  optimum  number 
greater  than  unity  of  dwells  per  scan  when  the  energy  allocated  for  target 
detection  is  fixed  (i.e.,  constant  m  *  n)  and  when  the  radar  transmission 
frequency  is  changed  from  dwell -to-dwel 1 .  Although  these  results  were 
first  derived  by  Swerling  (1960)  many  years  ago,  they  are  cast  here  in  a 
form  that  includes  both  coherent  and  noncoherent  integration. 

Then  three  cases  are  examined.  First,  the  case  that  the  propa¬ 
gation  channel  bandwidth  is  very  large  with  respect  to  the  radar  hopping 
bandwidth  (u>hop)  is  considered.  In  this  case  the  received  signal  con¬ 
tribution  from  the  fading  channel  remains  constant  during  the  entire  scan. 
Second,  results  are  presented  for  the  case  that  the  propagation  channel 
bandwidth  is  small  with  respect  to  the  radar  hopping  bandwidth.  In  this 
case  the  contribution  to  the  received  signal  from  scintillation  effects  is 
independent  from  dwel 1 -to-dwel 1 .  Third,  the  intermediate  situation  where 
the  channel  fluctuations  are,  correlated  from  dwell -to-dwel 1  is  considered. 
In  all  three  cases  here,  the  SBR  frequency  is  changing  between  dwells,  so 
that  a  Swerling  2  target  model  is  appropriate. 

In  the  following.  It  Is  shown  that  there  is  a  significant  gain 
in  SBR  detection  performance  for  the  case  that  the  propagation  scintilla¬ 
tion  effects  from  dwell -to-dwel!  are  totally  or  partially  decorrelated. 

Undisturbed  Environment 

Swerling  1  Target.  For  the  case  of  a  Swerling  1  target  and  no 
propagation  fading,  the  radar  cross  section,  and  therefore  the  target 
amplitude,  is  constant  for  all  the  dwells  of  a  scan.  For  a  fixed  target 
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amplitude,  the  probability  density  function  of  the  output  of  the  coherent 
integrator  is  given  by  Equation  14.  The  probability  density  function  of 
x  (see  Figure  3)  is  then  the  Rician  power  distribution 

p(xiA)  =  -L-  exp[  -(x+n2A/2)/no2]  I0  (/2a7/<j^)  (21) 

noN 

where  A  *  a2  is  the  constant  target  power  per  pulse. 

The  distribution  of  the  output  amplitude  from  the  sum  of  m 
amplitudes  each  with  the  probability  density  function  of  Equation  21  is 
obtained  by  noting  that  the  probability  distribution  of  the  sum  of  inde¬ 
pendent  random  variables  has  as  its  characteri stic  function  the  product  of 
the  characteristic  functions  of  the  individual  random  variables.  The 
amplitudes  x4  (z  3  l,2,...,m)  are  independent  because  the  noise  compon¬ 
ents  of  the  amplitudes  are  independent.  The  target  power  A  represents  a 
constant  component  which  is  common, to  all  of  the  terms.  ' 

The  characteristic  function  is  the  Fourier  transform  of  p(x|A) 

so  that 

Ci  (q  [A )  ■/  e“9x  p  ( x  I A  )dx  (22) 

'o' 

where  q  is  an  imaginary  transform  variable  and  where  the  subscript  1 
refers  to  the  fact  that  Ci (q | A )  is  the  characteristic  function. of  one  of 
the  amplitudes  in  the  noncoherent  integration.  Performance  of  the 
indicated  integration  (Gradshteyn  and  Ryzhik,  1965)  gives 

exp{-n2Aq/[2(no2q+l)]} 

Cj,  (q | A )  *  - -  -  -  (23) 

no2q+l 

Since  the  amplitude  z  is  the  sum  of  m  independent  values  of  x,  the 
characteristic  function  of  p(z|A),  the  probability  density  function  of 
z  for  a  given  target  power,  is  then 
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Cm(q|A)  =  [C  i  (q | A )  ]m 


exp{-mn2Aq/[2(no2q+l)  j } 

=  - 5 - - -  (24) 

(najjjq+1  )m 

The  probability  density  function  of  z  for  a- given  value  of  A,  p(z|A), 
is  the  inverse  Fourier  transform  of  Cm(q|A).  This  can  then  be  averaged 
over  the  scan-to-scan  distribution  of  A  to  obtain  p(z)  for  a  Sterling  1 
target.  Noting  that  the  order  of  integration  is  unimportant,  it  is  con¬ 
venient  to  average  Cm(q|A)  over  the  distribution  of  A  and  then  take 
the  inverse  Fourier  transform  to  obtain  p(z).  For  a  Swerling  1  target, 
the  probability  density  function  of  the  received  power  is 

p (A)  «  exp(-A/S0)/S0  (25) 

where,  again.  So  is  the  mean  received  power  per  pulse.  The  average  of 
Cm(q | A )  gives 


cm(q)  3  /  Cm(q|A)p(A)dA 


(no2q+l )ffl-l{l+no2q[l+mnS0/(2o^)  ]} 


The  inverse  Fourier  transform  (Campbell  and  Foster,  1948)  of  Cm(q)  gives 
the  probability  density  function  of  z  as  a  function  of  the  mean  signal- 
to-noise  ratio  as 

f  l+m<SNR>1m"1  exp{ -z/ [  no*  (l+m<SNR>)]f 


m<SNR> 


no2(l+m<SNR>) 


,  y{m-l,m<SNR>z/[  nq^(  l+m<SNR )] } 
r(m-l) 


where  y  is  the  incomplete  gamma  function,  and  r(m-l)  is  the  gamma 
function  (Abramowitz  and  Stegun,  1965).  For  later  reference  note  that 
r(a,x)  +  r(a,x)  *  r(a)  where  r(a,x)  is  also  referred  to  as  the  incom¬ 
plete  gamm.  function.  The  factor  <SNR>  is  the  mean  signal-to-noise 
ratio  per  dwe  1  defined  in  connection  with  Equation  19. 

Swerling  2  Target.  In  the  case  of  a  Swerling  2  target,  the 
radar  cross  section  is  independent  from  dwell -to -dwe 11,  and  the  proba¬ 
bility  density  function  of  the  amplitude  w  at  the  output  of  the  coherent 
integrator  is  given  by,  Equation  17  for  a  constant  value  of  S,  the  con¬ 
tribution  of  propagation  fading  to  the  received  power.  Note  that  Equation 
17  already  includes  the  contribution  of  the  Swerling  2  target. 

The  probability  density  function  of  x  then  has  the  exponential  form 

exp{-x/ [nojjj(l+S<SNR>)  ] }  1 

p(x|S)  » - 5—2 -  (28) 

1  nojj  (1+S<SNR>) 

where  <SNR>  is  the  mean  signal-to-noise  ratio  per  dwell.  Noting  that  the 
sum  of  m  independent  exponentially  distributed  random  variables  is  a  x2 
variable  with  2m  degrees  of  freedom,  the  probability  density  function  of 
z  for  a  fixed  value  of  S,  the  contribution  from  fading,  is 

iz/[noj£(l+S<SNR>)  }m-i  exp{-z/[no?(l+S<SNR>) ]} 

p(z|S) - 2 - - - : - - - - -  (29) 

1  no^(l+S<SNR>) T(m) 

where  T(m)  *  (m-1 ) !  Is  the  gamma  function. 

Results.  After  the  noncoherent  process  is  completed,  the 
amplitude  z  is  compared  with  a  threshold  and  a  hit  or  miss  is  declared. 
The  threshold  Is  set  so  that  the  probability  of  a  hit  from  noise  alone  is 
small  (typically  10”4  to  10"8).  The  probability  of  false  alarm  Is  then 
just  the  Integral  of  Equation  29  with  the  target  cross  section  set  to 


zero. 


=  /  p(zjS=0)dz  =  r[mJ/(noJ)]/r(m) 


where  r(m,x)  is  the  incomplete  gamma  function.  It  is  an  easy  matter  to 
numerically  invert  Equation  30  to  find  T/na^  as  a  function  of  the 
probability  of  false  alarm  and  the  number  of  dwells  per  scan,  m.  Since 
the  noise  is  independent  of  target  or  channel  statistics.  Equation  30  for 
the  threshold  does  not  depend  on  the  target  or  channel  model. 

Once  the  threshold  has  been  obtained,  the  probability  of  detec¬ 
tion  becomes 

Pd(S)  *  /  p(z|S)dz  (31) 

T 

which  is  the  probability  that  an  amplitude  z  of  signal  plus  noise 
exceeds  the  threshold  for  a  given  value  of  S,  the  signal  contribution  due 
to  fading. 

For  the  case  of  a  Swerling  1  target  with  no  fading,  it  is  neces- 
tary  to  perform  the  integration  given  by  Equation  31  using  Equation  27  as 
the  Integrand.  After  some  manipulation  the  probability  of  detection  is 
found  as 

pd  ■  pfa  *  exp(-T/[nojJ(lH«<SNR>)  ]} 


*  T{m-l,m<SNR>T/[no^(l+fn<SNR>)  ]} 

r(m-l) 
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Since  T/no^  is  independent  of  n,  the  probability  of  detection  depends 
only  on  the  mean  signal-to-noise  ratio  per  dwell,  the  number  of  dwells  per 
scan,  and  the  probability  of  false  alarm.  Also,  it  is  easy  to  show  that 
pd  =  pfa  when  <SNR>  is  zero. 

The  probability  of  detection  for  a  Swerling  2  target  has  the 
simpler  form 

pd(S)  *  r{m,T/[noJ(H*S<SNR>)]}/r(m)  (33) 

which  is  again  independent  of  n  for  a  given  signal -to-noise  ratio  per 
dwell.  Equation  33  gives  the  probability  of  detection  of  a  Swerling  2 
target  for  a  fixed  propagation  scintillation  condition  or  a  fixed  value  of 
S.  The  result  for  the  case  of  no  fading  is  found  by  taking  S  as  unity 
in  the  equation.  The  results  for  the  case  of  propagation  scintillation 
may,  in  some  cases  as  follows,  be  obtained  from  Equation  33  after  aver¬ 
aging  over  the  probability  distribution  function  of  S. 

The  probability  of  detection  is  shown  in  Figure  6  for  a  Swerling 
1  target  versus  the  mean  signal-to-noise  ratio  per  scan,.  m<SNR>,  The 
number  of  dwells  noncoherent ly  integrated  is  1,  2,  4,  8,  or  16.  In  the 
absence  of  target  cross  section  fluctuations  from  dwell -to-dwell ,  the 
signal-to-noise  ratio  required  to  achieve  a  given  probability  of  detection 
increases  as  the  number  of  dwells  per  scan  is  increased.  Thus  these  re¬ 
sults  demonstrate  the  well  known  result  that  for  a  Swerling  1  target  or, 
equivalently,  for  a  radar  that  does  not  change  transmission  frequencies 
from  dwell -to-dwell ,  the  optimum  number  of  dwells  per  scan  is  one. 

However,  when  the  radar  transmission  frequency  Is  changed  from 
dwell -to-dwell  so  that  the  target,  cross  section  fluctuations  are  described 
by  the  Swerling  2  model,  the  detection  corves  of  Figure  7  are  obtained. 


MEAN  SNR  (dJB/SCAN) 


Figure  6.  Probability  of  detecting  a  Sterling  1  target  in  an  undisturbed 
propagation  environment. 


Figure  7.  Probability  of  detecting  a  Swerling  2  target  In  an  undisturbed 
propagation  environment. 


,32 


From  the  figure,  for  signal -to-noi se  ratios  less  than  about  15  dB  per 
scan,  the  number  of  dwells  per  scan  which  minimizes  the  signal-to-noise 
ratio  required  to  achieve  a  given  probability  of  detection  depends  on  the 
probability  of  detection.  If  an  SBR  could  always  operate  at  signal-to- 
noise  ratios  above  15  dB  per  scan,  then  the  optimum  number  of  dwells  per 
scan  is  the  largest  number  possible.  In  general,  however,  it  is  desirable 
that  the  detection  performance  of  the  radar  degrade  as'  gracefully  as 
possible  as  the  signal-to-noise  ratio  falls  until  the  point  where  the 
probability  of  detection  falls  below  the  minim  in  (say  50  percent)  required 
to  maintain  a  radar  track.  At  a  50  percent  probability  of  detection,  the 
two  dwell  per  scan  case  requires  a  signal-to-noise  ratio  of  only  12.5  dB 
per  scan  whereas  the  16  dwell  per  scan  case  requires  14.5  dB  per  scan. 

The  optimum  number  of  dwells  per  scan  will  therefore  be  arbitrarily 
defined  by  finding  the  detection  curve  which,  for  probabilities  of 
detection  above  0.5,  minimizes  the  difference  between  the  signal-to-noise 
ratio  that  it  requires  to  achieve  a  given  probability  of  detection  and  the 
minimum  signal-to-noise  ratio  required  to  achieve  that  same  probability  of 
detection.  For  the  detection  curves  plotted  in  Figure  7,  the  optimum 
number  of  dwells  per  scan  is  4, 

By  comparing  Figures  6  and  7,  it  can  be  seen  that  for  probabili¬ 
ties  of  detection  greater  than  about  0.3  and  for  the  number  of  dwells  per 
scan  greater  than  unity,  a  Swerllng  2  target  requires  a  smaller  sfgnal-to- 
noise  ratio  to  achieve  a  given  probability  of  detection  than  does  a 
Swerling  1  target. 

The  results  presented  above  show  that  the  optimal  number  of 
dwells  per  scan  Is  about  4  when  the  transmission  frequency  Is  changed  from 
dwel 1-to-dwei 1 .  Otherwise,  a  single  dwell  per  scan  is. optimal.  However, 
these  waveform  design  considerations  are  limited  by  many  practical  con¬ 
straints.  The  maximum. duration  of  a  coherent  dwell  Is  limited,  for  exam¬ 
ple,  by  the  time  duration  of  radar  cross  section  fluctuations  caused,  by 


target  motion  or  by  other  considerations  such  as  the  coherence  time  of  the 
receiver  phase  reference.  On  the  other  hand,  the  minimum  number  of  pulses 
per  dwell  or  equivalently  the  minimum  duration  of  a  dwell  is  limited,  for 
example,  by  doppler  resolution  and  clutter  attenuation  requirements. 
Therefore,  it  is  not  possible  to  choose  the  number  of  dwells  per  scan 
merely  on  the  basis  of  target  detectability  considerations.  In  addition, 
an  SBR  may  be  required  to  operate  at  lower  probabilities  of  detection  than 
ground  based  radars  because  of  power  constraints.  Decreasing  probability 
of  detection  with  decreasing  signal-t:o-noise  ratio  might  then  drive  the 
waveform  design  to  2  or  3  dwells  per  scan. 

Effects  of  Fading  on  Noncoherent  Integration 

Now  consider  the  case  that  the  signal  decorrelation  time  is 
large  with  respect  to  the  duration  of  all  the  dwells  that  comprise  a  scan. 
Under  this  slow  fading  assumption,  the  propagation  effects  are  controlled 
by  the  relationship  between  the  SBR  hopping  bandwidth,  uhop»  and  the 
prooagation  channel  coherence  bandwidth,  wCOh.  In  the  following,  the 
radar  frequency  Is  changed  for  each  dwell  So  that  a  Swerling  2  target 
model  applies. 

Large  Channel  Coherence  Bandwidth.  For  the  case  that 
(*>coh>  ^hop ’  contribution  to  the  received  signal  due  to  channel 
fluctuations  is  constant  during  a  scan  but  varies  independently  from  scari- 
to-scan.  Thus,  the  probability  of  detection  during  fading  is  computed  by 
averaging  the  probability  of  detection  for  constant  propagation  conditions 
over  the  appropriate  distribution  function  describing  the  propagation 
channel  fluctuations.  That  Is, 

vfpd(s)v°i>(s)  as  ■  1341  . 
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where  pprc,p(S)  is  the  probability  density  function  for  one-way,  mono- 
static,  or  bistatic  propagation  geometries.  In  general,  it  is  necessary  . 
to  evaluate  Equation  34  using  numerical  quadrature  techniques. 


Small  Channel  Coherence  Bandwidth.  If  u 


"“hop*  the  contn  ~ 


bution  to  the  signal  from  propagation  scintillation  .is  independent  from 
lwel 1-to-dwel 1 .  In  this  case,  the  distribution  of  the  amplitude  out  of 
the  square-law  detector  (Equation  33)  must  first  be  averaged  over  the 
dwel 1 -to-dwel 1  distribution  of  S  before  the  noncoherent  integration  pro¬ 
cess  can  be  considered.  To  the  authors'  knowledge,  this  average  cannot  be 
obtained  in  closed  form  for  any  of  the  propagation  geometries  considered 
here.  It  is  therefore  nec^su.  ■  to  use  Monte  Carlo  techniques  to  obtain 
the  probability  of  detection. 


Here  random  samples  of  S  are  generated  from  the  appropriate 
distribution  and  then  Equation  28  is  used  to  generate  m  random  samples 
of  the  output  amplitude  from  the  square-law  detector.  The  amplitudes  are 
noncoherent ly  integrated  and  compared  with  the  threshold  from  Equation  30 
to  declare  a  hit  or  a  miss.  The  probability  of  detection  is  then  esti¬ 
mated  as  the  ratio  of  the  number  of  hits  to  the  number  of  thials  or  scans. 

The  samples  of  S,  the  fading  propagation  channel  contribution 
to  the' signal,  are  generated  by  noting  that  for  one-way  propagation  the 
probability  density  function  of  S  Is  exponential: 

P(S)  -  e-S  .  (35) 

Independent  samples  of  S  are  then  obtained  by  Inverting  the  equation 


f{$)  -  /  p(S '  )dS'  »  1  -  e's  ■  i-U 

o 


(36) 


where  F(S)  is  the  cumulative  distribution  of  S  and  U  is  a  uniformly 
distributed  random  variable  on  the  interval  (0,1].  For  a  monostatic  SBR 
geometry,  S  is  the  square  of  the  one-way  value;  for  the  bistatic  SBR 
geometry,  S  is  the  product  of  two  independent  one-way  values'.  Hence  the 
values  of  S  are  generated  in  all  three  SBR  geometries  as 

-In  U  one-way 

i.  in2U  monostatic  (37) 

2 

in  U  in  U1  bistatic 

where  U  and  U:  are  independent  variates  with  uniform  distributions. 

The  mean  value  of  S  is  unity  in  each  case. 

For  a  given  value  of  S  generated  from  Equation  37,  the  samples 
of  the  amplitude  at  the  output  of  the  square-law  detector  are  obtained  by 
inverting  the  cumulative  distribution  of  amplitude.  In  this  case  the 
cumulative  distribution  function  of  the  output  amplitude  x  is  given  as 
the  integral  of  Equation  28  with  respect  to  x.  Random  samples  of  the 
amplitude  may  then  be  obtained  in  the  same  manner  as  above  from  the 
expression 

x/noj  =  -(1  +  <SNR>S)tn  U"  (38) 

where  <SNR>  is  the  input  value  of  the  mean  signal-to-noise  ratio  per 
dwell  and  where  U’1  is  another  uniformly  distributed  random  variable. 

The  noncoherently  integrated  amplitude  may  then,  be  written  as 
m 

z/no2  *  -  l  (l  +  <SNR>S.)fcn  U. 

,  N  1-1  .  1  1 


(39) 


where  independent  values  of  S  and  U  are  generated  for  each  term  of  the 
sum.  The  amplitude  z/nu^  is  then  compared  with  the  threshold  T/naj*  to 
declare  a  hit  or  a  miss.  Note  that  both  z/noj*  and  T/noj^  are  indepen¬ 
dent  of  the  values  of  n  and  so  that  the  probability  of  detection  is 

only  a  function  of  the  probability  of  false  alarm,  the  mean  signal-to- 
noise  ratio  per  dwell,  and  the  number  of  dwells  per  scan. 

The  accuracy  of  the  Monte  Carlo  technique  can  be  checked  by 
generating  the  probability  of  detection  for  the  case  of  a  single  dwell  per 
s,can  and  comparing  the  result  with  that  calculated  from  Equation  34.  In 
this  case  the  results  for  both  small  and  large  coherence  bandwidth  should 
be  identical. 

General  Channel  Coherence  Bandwidth.  Intermediate  between  the 
large  and  the  small  channel  coherence  bandwidth  cases  are  propagation 
environments  where  the  dwel 1 -to-dwel 1  transmission  frequency  changes  par¬ 
tially  decorrelate  the  propagation  effects.  Slow  fading  is  again  assumed 
so  that  the  propagation  effects  are  constant  for  the  duration  of  a  dwell. 
Once  the  m  correlated  samples  S-j  have  been  obtained.  Equation  39  is 
used  to  calculate  the  noncoherently  integrated  amplitude  z/no^.  it  is 
convenient  in  this  situation  to  generate  the  correlated  samples  of  S  for 
the  three  SBR  geometries  from  the  underlying  normally  distributed  one-way 
voltage. 

It  is  known  (Fante,  1975)  that  after  propagation  one-way  through 
a  strongly  turbulent  layer,  the  received  voltage  of  an  initially  constant 
amplitude  and  phase  signal  can  be  expressed  as  the  sum  of  quadrature 
components 

V*  V  <«» 
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where  is  the  complex  received  voltage  on  the  kth  dwell  with 
in-phase  component  I|<  and  quadrature-phase  component  Q^.  The 
components  of  Ej<  are  independent,  normally  distributed  variates  with 
zero  mean  and  variance  equal  to  1/2  so  that,  in  the  one-way  propagation 
case,  the  power 

sk,i  ■  ll  *  •  <«> 

has  an  exponential  probability  density  function  with  mean  of  unity.  . 

A  simple  method  to  obtain  m  correlated  samples  of  E|<  is  to 
generate  a  Markov  process  (Mitchell,  1976) 

•  El’-  Ci 

(42) 

Efc  »  p  Ek_j  +  /l-p2  ,  1  <  k  <  m 

where  { }  is  a  sequence  of  complex,  independent  samples  with  first- 
order  statistics  identical  to  those  of  {E|<}.  It  is  then  straight¬ 
forward  to  show  that  the  autocorrelation  function  of  the  sequence  (E^J 
is 

<EkE*>  -  plk-*l  .  (43) 

Once  the  correlated  sequence  {E^}  has  been  generated  with  statistics 
appropriate  to  the  one-way  propagation  case,  the  sequence  of  resulting 
values  of  power  is  then  available  for  the  monostatic  SBR  geometry.  An 
additional  sequence  is  necessary  for  the  bistatic  SBR  configuration. 

The  values  of  ,  1  <  k  <  m  ,  for  the  m  dwells  of  a  scan  are  then 
given  by 
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one-way 

monostatic  (44) 

bistatic 

where  ^  is  an  independent  power  appropriate  to  the  second  one-way 
path  of  the  bi static  geometry. 

Results.  Figures  8-10  show  the  probability  of  detection  during 
slow  fading  for  the  three  cases  of  <DCOh  under  consideration.  The 
results  are  shown  as  a  function  of  signal-to-noise  ratio  and  the  number  of 
dwells  per  scan,  m.  A  check  of  the  results  is  made  possible  by  comparison 
of  the  detection  curve  for  one  dwell  per  scan  to  the  appropriate  detection 
curves  for  slow  fading  shown  in  Figure  4.  As  can  be  seen,  the  results  are 
identical  in  all  cases.  Figures  8a-c  show  the  probability  of  detection 
for  one-way,  monostatic,  and  bistatic  propagation  geometries,  respective¬ 
ly.  Because  the  probability  of  detection  is  insensitive  to  the  number  of 
dwells  per  scan,  only  results  for  m  3  1,  4,  and  16  are  shown.  These 
results  and  those  for  2  or  8  dwells  per  scan  are  summarized  in  Table  1, 
which  gives  the  probability  of  detection  as  a  function  of  the  mean  signal- 
to-noise  ratio  per  scan.  It  can  be  seen  from  the  table  that  for  proba¬ 
bilities  of  detection  in  the  range  0.5  to  0.95  where  the  SBR  would  normal¬ 
ly  be  expected  to  operate,  the  optimum  detection  performance  is  achieved 
with  four  dwells  per  scan.  However,'  the  difference  is  slight  between  the 
cases  of  2  and  4  dwells  per  scan.  An  inspection  of  the  figures  shows 
that,  for  a  given  probability  of  detection  above  about  0.7,  the. one  dwell 
per  scan  waveform  requires  the  highest  signal-to-noise  ratio.  At  a 
probability  of  detection  of  0.9,  the  difference  in  the  required 
signal-to-noise  ratio  between  the  case  of  one  dwell  and  of  four  dwells  per 
scan  is  2.5  dB  for  the  one-way  and  bi static  geometries  and  is  2  dB  for  the 
monostatic  geometry. 
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Table  1.  Probability  of  detection*  under  slow  Rayleigh  fading  conditions 
with  a  large  frequency  coherence  bandwidth. 


a).  One-Way  Geometry 


Mean  SNR 

Dwells/Scan 

1 

2 

4 

8 

16 

dB/Scan 

0.00 

5.755E-03 

2.038E-03 

5.499E-04 

1. 189E-04 

2.453E-05 

5.00 

5.828E-02 

3.848E-02 

2. 102E-02 

9.139E-03 

3.031E-03 

10.00 

2.314E-01 

2. 172E-01 

1.838E-01 

1. 378E-01 

8.900E-02 

15.00 

4.916E-01 

5.245E-01 

5.204E-01 

4.835E-01 

4.203E-01 

20.00 

7.241E-01 

7.802E-01 

7.955E-01 

7.829E-01 

7.501E-01 

25.00 

8.698E-01 

9. 154E-01 

9.270E-01 

9.237E-01 

9. 117E-01 

30.00 

9.438E-01 

9. 704E-01 

9. 755E-01 

9. 749E-01 

9. 710E-01 

35.00 

9. 773E-01 

9.897E-01 

9.916E-01 

9.918E-01 

9. 907E-01 

40.00 

9.912E-01 

9.961E-01 

9.968E-01 

9.972E-01 

9.970E-01 

b),  Monostatic  Geometry 


c)  Bistatic  Geometry 
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Mpan  SNR 

Dwell  s/Scan 

dB/Scan 

1 

2 

4 

3 

16 

0.00 

1.389E-02 

9.401E-03 

5.730E-03 

3. 103E-03 

1.464E-03 

5.00 

6.137E-02 

5.299E-02 

4.218E-02 

3.075E-02 

2.032E-02 

10.00 

1.679E-01 

1.6516-01 

1.511E-01 

1.293E-01 

1.035E-01 

15.00 

3.264E-01 

3.395E-01 

3.311E-01 

3.068E-01 

2. 720E-01 

20.00 

5.017E-01 

5.294E-01 

5.287E-01 

5.091E-01 

4. 768E-01 

25.00 

6.585E-01 

6.914E-01 

6. 948E-01 

6.816E-01 

6.575E-01 

30.00 

7. 796E-01 

8.090E-01 

8. 131E-01 

8.051E-01 

7.892E-01 

35.00 

8.638E-01 

8.860E-01 

8.894E-01 

8.848E-01 

8. 751E-01 

40.00 

9. 184E-01 

9.334E-01 

9.357E-01 

9.333E-01 

9.276E-01 

Mean  SNR 

Dwells/Scan 

dB/Scdn 

.  1 

2 

4 

8 

16 

0.00 

1. 114E-02 

6.538E-03 

3. 363E-03 

1.495E-03 

5.620E-04 

5.00 

6.137E-02 

4. 957E-02 

3.640E-02 

2.405E-02. 

1.408E-02 

10.00 

1.897E-01 

1.831E-01 

1.635E-01 

1.350E-01 

1.027E-01 

15.00 

3.376E-01 

4.045E-01 

3.945E-01 

3.634E-01 

3. 180E-01 

20.00 

5.971E-01 

S.354E-01 

6.385E-01 

6. 176E-01 

5i 798E-01 

25.00 

7.657E-01 

8.078E-01 

8. 158E-01 

8.053E-01 

7.826E-01 

30.00 

8. 768E-01 

9.096E-01 

9. 162E-01 

9. 119E-01 

9.007E-01 

35.00 

9.402E-01 

9.607E-01 

9.646E-01 

9.631E-01 

9. 583E-01 

.40.00  . 

9. 727E-01 

9.836E-01 

— r? - 

9.855E-01 

9.853E-01 

9.835E-01 

9 


10 


30 


40 


lllllliilllllfcllii.i 
19  20  29 

MEAN  SNR  (dB/SCAN) 


■  ■  »  *  1  •  *  *  * 
JS 


Figure  8(a).  Probability  of  detecting  a  Sterling  2  target  for 
one-way  SBR  geometry, 


Figure  8(b). 


Probability  of  detecting  a  Sterling  2  target  for  a 
monostatic  SBR  geometry,  «*coh» 


41 


Figure  8(c).  Probability  of  detecting  a  Sterling  2  target  for  a. 
bi static.  S8R  geoaetry,  “coh> 

The  detection  curves  for  small  coherence  bandwidth  are  shown  in 
Figures  9a-c  for  one-way,  monostatic,  and  bistatic  propagation  geometries, 
respectively.  The  simulation  results  are  shown  as  dots  connected  by 
straight  lines.  For  one  dwell  per  scan,  close  agreement  is  seen  between 
analytical  results  (solid  line)  and  the  simulation  results  (stars).  These 
detection  curves  have  the  same  qualitative  appearance  as  those  for  a 
Swerling  2  target  in  the  absence  of  propagation  fading.  That  Is,  at 
signal-to-noise  ratios  per  scan  above  15  dB  or  so,  increasing  the  number 
of  dwells  per  scan  Increases  the  probability  of  detection  at  the  expense 
of  giving  a  less  graceful  degradation  of  the  probability  of  detection  as 
the  signal-to-noise  ratio  Is  reduced.  If  the  criterion  for  defining  the 
optimum  number  of  dwells  per  scan  that  was  used  earlier  for  a  Swerling  2 
target  in  an  ambient  environment  is  again  applied,  the  optimum  number  of 


Figure  9(a).  Probability  of  detecting  a  Sterling  2  target  for  a 

one-wy  9M»try,  -cot|« 


Figure  9(b).  Probability  of  detecting  a  Sterling  2  target  for  a 
aonostatlc  S8R  geometry,  «CQh<  ^  . 
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Figure  9(c).  Probability  of  detecting  a  Sterling  2  target  for  a 
bistatic  SB*  gewetry.  »cot,<  "hop- 

dwel  ls  per  scan  for  o>  .  <  to.  lies  between  4  and  8.  Of  course,  as  dis¬ 
con  nop 

cussed  previously,  other  design  considerations  might  force  the  use  of 
fewer  than  the  optimum  number  of  dwells  determined  on  the  basis  of 
detectability  alone. 

Figures  lOa-c  show  the  effects  of  partial  decorrelation  of  the 
propagation  channel  between  dwells  for  the  three  SBR  geometries.  The 
probability  of  detection  versus  the  mean  signal -to-noise  ratio  per  scan  is 
shown  for  8  dwells  per  scan  with  the  correlation  coefficient  defined  in 
Equation  42,  ranging  from  1  {which  corresponds  to  c«coh>  ‘“hop)  t0  ® 

(which  corresponds  to  <*>CQh«  ^op^'  Pr°bal)iHty  of  detection  with'  8 
dwells  per  scan  in  an  ambient  environment  is  also  shown.  These  results 
show  that  only  a  small  amount  of  inter-dwell  decorrelation  of  the  fading 
effects  is  necessary  to  substantially  improve  target  detectability  in 
comparison  to  the  case  of  no  decorrelation.  For  example,  at  a  probability 
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Figure  10(a).  Probability  of  detecting  a  Sterling  2  target  for  a  one-way 
SBR  geometry  characterized  by  dwel 1 -to-dwel 1  correlation. 


Figure  10(b).  Probability  of  detecting  a  Swerllng  2  target  for  a  mono 
static  S8R  geometry  characterized  by  dwell -to-dwel 1 
correlation. 
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Figure  10(c).  Probability  of  detecting  a  Sterling  2  target  for  a 

bi static  SBR  geometry  characterized  by  dMel 1 -to-dwel 1 
correlation. 

of  detection  of  0.7,  a  monostatic  SBR  that  utilizes  8  dwells  per  scan  will 
suffer  losses  in  detection  sensitivity  of  2.7  dB  for  p  *  0  ,  4.9  dB  for 
p  *  0.8  ,  and  11  dB  for  p  *  1.  While  it  is  apparent  that  the  use  of  mul¬ 
tiple  dwells  per  scan  does  not  improve  detection  performance  in  a  slow 
fading,  large  coherence  bandwidth  cha.mel,  only  a  small  amount  of  dwell- 
to-dwell  decorrelation  of  the  scintf+Vatton  effects  substantially  improves 
detection  performance.  ' 


SECTION  5 
CONCLUSIONS 


It  is  apparent  that  severe  fading  has  a  dramatic  effect  on  the 
target  detection  performance  of  An  SBR.  For  slow  *'ading,  when  there  is  no 
degradation  to  the  SBR  coherent  integration  process,  the  losses  in  detec¬ 
tion  sensitivity  are  7  dB  for  bistatic  operation  and  11  dB  for  monostatic 
operation.  For  fast  fading,  severe  degradation  occurs  when  the  signal  de- 
correlation  time  becomes  smaller  than  about  a  tenth  of  the  SBR  coherent 
integration  time.  A  design  alternative  useful  for  fast,  fading  might  in¬ 
volve  adaptively  changing  the  dwell  duration  so  that  the  SBR  coherent 
integration  time  does  not  exceed  the  signal  decorrelation  time. 

When  the  total  signal-to-noise  ratio  per  scan  for  the  detection 
of  targets  is  fixed,  the  optimum  detection  performance  in  ambient  or  slow 
fading  conditions  is  achieved  when  the  radar  energy  is  divided  into  2  to  8 
dwells  per  scan.  Each  dwell  rjiust  use  a  separate  transmission  frequency  so 
that  Swerling  2  radar  cross  section  statistics  apply.  Utilization  of  more 
dwells  per  scan  enhances  the  detection  performance  at  high  signal-to-noise 
ratios  at  the  expense  of  detection  curves  that  degrade  less  gracefully  as 
the  signal  power  decreases. 

When  multiple  dwell;  per  scan  are  used  during  slow  fading,  the 
detection  performance  is  sensitive  to  the  coherence  bandwidth  of  the  prop¬ 
agation  channel.  For  example,  for  an  8  dwell  per  scan  waveform  and  at  a 
probability  of  detection  of  0,7,  the  difference  in  detection  sensitivity 
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between  the  worst  case  large  coherence  bandwidth  environment  where  propa¬ 
gation  effects  are  constant  during  the  scan,  and  the  best  case  small 
coherence  bandwidth  environment  where  scinti 1 lation,  effects  give 
independent  voltages  from  dwell -to-dwell ,  is  8  dB  for  a  monostatic  radar 
and  5.5  dB  for  a  bistatic  radar.  However,  only  a  small  amount  of 
dwell -to-dwell  correlation  of  the  scintillation  is  necessary  to  achieve 
most  of  this  detection  sensitivity  gain.  It  is  then  apparent  that  an  SBR 
should  attempt  to  maintain  a  large  enough  transmission  frequency 
separation  between  dwells  to  insure  that,  the  scintillation  effects 
partially  decorrelate  from  dwell -to-dwell .  Hence,  an  SBR  measurement  of 
the  channel  bandwidth  appears  a  useful  mitigation  tool. 

Both  mitigation  techniques  suggested  here  will  affect  other 
aspects  of  SBR  performance  and  require  further  investigation  in  the 
context  of  a  complete  SBR  system. 
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